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ABSTRACT 



> 

Based on Cosmic Microwave Background (CMB) maps from the 2013 Planck Mission data release, this paper presents the detection of the 
ly-^ Integrated Sachs-Wolfe (ISW) effect, i.e., the correlation between the CMB and large-scale evolving gravitational potentials. The significance of 
• detection ranges from 2 to 4 cr, depending on which method is used. We investigate three separate approaches, which cover essentially all previous 
studies, as well as breaking new ground, (i) Through correlation of the CMB with the Planck reconstructed gravitational lensing potential (for the 
first time). This detection is made using the lensing-induced bispectrum between the low-f and high-f temperature anisotropics; the correlation 
between lensing and the ISW effect has a significance close to 2.5 cr. (ii) Through cross-correlation with tracers of large-scale structure, yielding 
around 3 cr significance, based on a combination of radio (NVSS) and optical (SDSS) data, (iii) Using aperture photometry on stacked CMB fields 
at the locations of known large-scale structures, which yields and confirms, over a broader spectral range, a 4 cr signal when using a previously 
explored catalogue, but shows strong discrepancies in amplitude and scale compared to expectations. More recent catalogues give more moderate 
results, ranging from negligible to 2.5 cr at most, but with a more consistent scale and amplitude, the latter being still slightly above what is 
< I expected from numerical simulations within ACMD. Where they can be compared, these measurements are compatible with previous work using 
data from WMAP, which had already mapped these scales to the limits of cosmic variance. Planck's broader frequency coverage allows for better 
foreground cleaning, and confirms that the signal is achromatic, bolstering the case for ISW detection. As a final step we use tracers of large-scale 
structure to filter the CMB data, presenting maps of the ISW temperature perturbation. These results provide complementary and independent 
evidence for the existence of a dark energy component that governs the current accelerated expansion of the Universe. 

Key words. Cosmology: observations ~ cosmic microwave background - large-scale structure of the Universe - dark engery - Galaxies: clusters: 
general - Methods: data analysis 
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1. Introduction 

This paper, one of a set associated with the 2013 data re- 
lease from the Planck^ mission (Planck Collaboration I 2013), 
presents the first results on the integrated Sachs-Wolfe (ISW) ef- 
fect using Planck data. The ISW effect (Sachs & Wolfe 1967; 
Rees & Sciama 1968; Hu & Sugiyama 1994) is a secondary 
anisotropy in the cosmic microwave background (CMB), caused 
by the interaction of CMB photons with the time-evolving po- 
tentials from large-scale structure (LSS, hereafter). Photons fol- 
low a geodesic that is weakly perturbed by the Newtonian grav- 
itational potential, O, and experience a fractional shift in their 
temperature given by 




where the integral is expressed in terms of the conformal time 
77, defined differentially by drjlda - l/(a^H(a)) with H{a) the 
Hubble function and a the scale factor. The integration limits 
here go from the recombination time (77*) to the present time 

(m)- 

The sensitivity of the ISW effect to gravitational potentials 
(that can extend over Gpc scales) results in the power of the 
ISW being concentrated on the largest scales. The largest scales 
for the CMB have been mapped out by the Wilkinson Microwave 
Anisotropy Probe (WMAP) to the statistical limit of cosmic vari- 
ance. Some systematics (like foreground removal) can have an 
impact on the reconstruction of the CMB especially at the largest 
scales where our Galaxy can introduce significant residuals on 
the reconstructed CMB map. The superior sensitivity of Planck 
together with its better angular resolution and wider frequency 
coverage allows for a better understanding (and hence removal) 
of Galactic and extagalactic foregrounds, therefore reducing the 
possible negative impact of these residuals. Planck allows us to 
improve on previous measurements by having a better system- 
atic control, an improved removal of foregrounds (that permits 
us to explore the achromatic nature of the ISW signal on a wider 
frequency range), and a better understanding of systematics af- 
fecting tracer catalogues. 

For cosmological models where Q.^ = 1, gravitational po- 
tentials remain constant during linear structure formation, and 
the ISW signal is negligible (to first order, although second or- 
der nonlinear ISW is always expected around smaller over- and 
under-dense regions). In the presence of dark energy, decaying 
potentials due to the accelerated expansion rate, result in a net 
ISW effect which is positive when the CMB photons cross over- 
dense regions and negative when the CMB photons cross under- 
dense regions. Therefore, the ISW effect is an indicator of either 
non-zero curvature (Kamionkowski & Spergel 1994), any form 
of dark energy, such as a cosmological constant A (Crittenden & 
Turok 1996), modified gravity (Hu 2002), or a combination of 
these possibilities. By measuring the rate at which gravitational 
potentials in the LSS decay (up to redshift of around 2), the ISW 
effect can be used as an independent probe of cosmology and 
provides complementary and independent evidence for dark en- 
ergy. 

Detection of the ISW effect was first made possible with 
all-sky CMB maps from WMAP. Based on these data, many 

' Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



works can be found in the literature where the authors aim at 
making, and subsequectly improving, the measurement of the 
ISW effect through correlations with tracer catalogues: 2MASS 
(an infrared catalogue at low redshifts around 0.1, Afshordi 
et al. 2004; Dupe et al. 2011; Francis & Peacock 2009; Rassat 
et al. 2006), HEAO (an X-ray survey at low redshift, with the 
first positive claim for detection, Boughn & Crittenden 2004), 
Sloan Digital Sky Survey (SDSS, optical survey at intermediate 
redshifts, Bielby et al. 2010; Cabre et al. 2006; Fosalba et al. 
2003; Fosalba & Gaztaiiaga 2004; Giannantonio et al. 2006; 
Granett et al. 2009; Lopez -Corredoira et al. 2010; Padmanabhan 
et al. 2005; Sawangwit et al. 2010; Scranton et al. 2003; Xia 
2009), the NRAO VLA Sky Survey (NVSS, radio catalogue with 
high-redshift sources, Boughn & Crittenden 2005; Hernandez- 
Monteagudo 2010; Massardi et al. 2010; McEwen et al. 2007; 
Pietrobon et al. 2006a; Raccanelli et al. 2008; Schiavon et al. 
2012; Vielva et al. 2006), and combined measurements with 
multiple tracers (Corasaniti et al. 2005; Gaztanaga et al. 2006; 
Giannantonio et al. 2008, 2012; Ho et al. 2008; Nolta et al. 
2004). The significance of the ISW detections that can be found 
in the literature range between 0.9 cr and 4.7 cr. There are a num- 
ber of peculiarities related to some of the detection claims, as 
noted by Hernandez -Monteagudo (2010) and Lopez -Corredoira 
et al. (2010). They both found lower significance levels than 
some previous studies and pointed out the absence of the signal 
at low multipoles where the ISW effect should be most promi- 
nent and the presence of point source emission on small scales 
for radio surveys. 

The main result that is obtained from an ISW detection is 
a constraint on the cosmological constant (or dark energy), Oa. 
The general consensus from the variety of ISW analyses is for 
a value of Qa ^ 0.75 with an error of about 20%, which pro- 
vides independent evidence for the existence of dark energy 
(Fosalba et al. 2003; Fosalba & Gaztanaga 2004; Nolta et al. 
2004; Corasaniti et al. 2005; Padmanabhan et al. 2005; Cabre 
et al. 2006; Giannantonio et al. 2006; Pietrobon et al. 2006b; 
Rassat et al. 2006; Vielva et al. 2006; McEwen et al. 2007; Ho 
et al. 2008). All tests on spatial flatness find an upper limit for 
D.K of a few percent (Nolta et al. 2004; Gaztaiiaga et al. 2006; 
Ho et al. 2008; Li & Xia 2010). Using a prior on spatial flat- 
ness, the dark energy equation of state parameter, w, was found 
to be close to -1 (Giannantonio et al. 2006; Vielva et al. 2006; 
Ho et al. 2008) and has been excluded from having a strong time 
evolution (Giannantonio et al. 2008; Li & Xia 2010). 

The ISW effect is achromatic, conserving the Planck spec- 
trum of the CMB and can be separated from other CMB fluc- 
tuations through cross-correlations with catalogues which trace 
the LSS gravitational potentials (see for instance Crittenden & 
Turok 1996). This cross-correlation can be studied in different 
ways: angular cross-correlations in real space between the CMB 
and the catalogues tracing the LSS; the corresponding angu- 
lar cross-power spectrum of the Fourier-transformed maps; or 
through the covariance of wavelet-filtered maps as a function of 
wavelet scale. The studies using WMAP data mentioned above 
follow this survey cross-correlation techique. 

An alternative approach, similar to the angular cross- 
correlation in real space, consists of stacking CMB fields centred 
on known supersclusters or supervoids (Granett et al. 2008a,b; 
Papal & Szapudi 2010). The advantage of this technique is that 
it allows for a detailed study of the profile of the CMB fluctua- 
tions caused by this secondary anisotropy. 

A novel and powerful approach takes advantage of the fact 
that the same potentials that make CMB photons gain or loose 
energy along their path (ISW), create lensing distortions that can 
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be measured from the CMB map directly (e.g., Hu & Okamoto 
2002). The interplay between weak gravitational lensing and 
the ISW effect causes a non-Gaussian contribution to the CMB, 
which can be measured through the lensing-induced bispectrum 
between small and large angular scales. The measurement of the 
lensing potential requires a large number of modes that could not 
be measured before the arrival of Planck data. 

This paper presents new measurements of the ISW effect car- 
ried out with Planck. Even although our detections are not in ev- 
ery case as strong as some previously claimed significance lev- 
els, we believe that our results are an improvement over earlier 
studies. This is because we can use the additional power enabled 
by the frequency coverage and sensitivity of Planck. To estab- 
lish this we carry out a comprehensive study of all the main 
approaches which have previously been taken to estimate the 
ISW signal. We also present new results in relation to the non- 
Gaussian structure induced by the ISW effect. 

The paper is organized as follows: In Sect. 2 we describe 
the data used in this work (both for the CMB and large-scale 
structure). The first ever results on the estimation of the lensing- 
induced bispectrum are presented in Sect. 3. Cross-correlations 
with external surveys are investigated in Sect. 4, and in Sect. 5 
we present the results for the stacking analysis on the tempera- 
ture maps, as well as aperture photometry on super-clusters and 
super-voids. The recovery of the ISW all-sky map is described 
in Sect. 6. Finally, we discuss our main results and their cosmo- 
logical implications in Sect. 7. 

2. Data description 

In this Section we describe the different data sets used in this 
paper. This includes Planck data (the CMB temperature and 
lensing potential maps, see Planck Collaboration I 2013; Planck 
Collaboration XII 2013; Planck Collaboration XVII 2013) we 
refer to the corresponding Planck papers for details) and exter- 
nal data sets (large-scale structure tracers) used in the ISW deter- 
mination: the radio NVSS catalogue; optical Luminous Galaxies 
(CMASS/LOWZ) and the Main Galaxy Sample from the Sloan 
Digital Sky Survey (SDSS); as well as several superstructure cat- 
alogues. 

2.1. Planck data 

Planck data and products used in this paper are described in the 
following sections, in particular the foreground-cleaned CMB 
maps produced by the Planck component separation pipelines, 
and related products, such as dedicated component-separated 
frequency maps (Planck Collaboration XII 2013), as well as the 
Planck lensing map (Planck Collaboration XVII 2013). 

2.1.1. CMB maps 

For the present work we have made use of the Planck 
foreground-cleaned CMB maps provided by the data process- 
ing centres (as described in the Planck component separation 
paper Planck Collaboration XII 2013). In particular, to test ro- 
bustness, some of the results are presented for different cleaned 
CMB maps, which were constructed using four different com- 
ponent separation techniques: Coraraader-Ruler (C-R, which 
uses physical parametrization); NILC (an internal linear combi- 
nation technique); SEVEM (a template fitting method); and SMICA 
(which uses spectral matching). Since the contribution of the 
ISW signal is only important on large scales, low resolution 



maps, with HEALPix Gorski et al. (2005) parameter A^side = 64, 
and pixel size of about 55 arcmin, have been used for most of the 
analyses. One exception is the study of the correlation between 
the ISW and lensing signals, which requires the use of full- 
resolution maps (A^side=2048, pixel size of 1.7 arcmin). The maps 
are degraded directly from the original full resolution down to 
the corresponding A^side- 

In addition, foreground-cleaned maps per frequency (from 
44 to 353 GHz) at resolution A^side = 512 were used for the stack- 
ing analysis presented in Sect. 5. These maps were constructed 
by subtracting a linear combination of internal templates using 
SEVEM (see the SEVEM Appendix of Planck Collaboration XII 
2013, for a detailed description of the method). As an example 
the SEVEM CMB map is shown in Fig. 1 (left panel). 

Finally, to minimize the presence of foreground contami- 
nation in the maps, we have used the official mask described 
in Planck Collaboration XII (2013), which excludes regions with 
larger Galactic and point-source contamination (the U73 mask). 
This mask is given at the full Planck resolution and is down- 
graded to the required levels. The downgrading procedure con- 
sists of the following steps: the mask (originally a map with zero 
and one values) is convolved with a Gaussian beam of FWHM 
three times the characteristic pixel size of the final A^side resolu- 
tion; this convolved map is then degraded to the required A^side. 
and, afterwards, a threshold of 0.75 is imposed (i.e., pixels with 
a value above this threshold are set to one, whereas the rest are 
set to zero). 

2.1.2. Lensing potential map 

Weak gravitational lensing distorts the CMB temperature 
anisotropy pattern. This effect is sensitive to the projected matter 
distribution in the large-scale structure at high redshifts, where 
structure growth is linear and the statistics close to Gaussian. 
Weak lensing causes correlations between different multipoles 
which are proportional to the lensing deflection field. These cor- 
relations can be exploited for reconstructing the density field and 
for measuring its statistical properties (Hu & Okamoto 2002; 
Okamoto & Hu 2003). The lensing effect in the CMB can be 
estimated by this homogeneity breaking, and in this way indi- 
vidual modes of the lensing potential at multipoles f < 100 can 
be reconstructed with a significance of around 0.5 cr, showing 
the necessity of a statistical treatment. Nevertheless the over- 
all effect of the lensing is measured to better than 25 cr (Planck 
Collaboration XVII 2013). The additional lensing effect in the 
temperature power spectrum is detectable with a significance of 
about 10 o- (Planck Collaboration XV 2013). 

With Planck data, we aim at detecting a correlation between 
the ISW effect and the lensing potential, where the latter is a 
tracer of the large-scale structure at high redshift. This correla- 
tion is restricted to 9 cr, even in the ideal case, limited by cos- 
mic variance and the smallness of the ISW effect in compari- 
son to the primary CMB (Lewis et al. 2011). The data products 
used in this study are the Planck lensing potential reconstruction, 
and specific lensing maps obtained from the component separa- 
tion pipelines. The lensing potential is available as part of the 
first Planck data release. Its detailed development is described in 
the Planck lensing paper (Planck Collaboration XVII 2013). In 
Fig. 1 we reproduce (right panel) an optimally filtered version 
of the Planck lensing map, suitable for the ISW-lensing cross- 
correlation. 

In addition to a direct correlation between the CMB sky and 
the reconstructed lensing map, we measure the bispectrum gen- 
erated by weak lensing by applying a range of estimators: the 
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Fig. 1. Left: one of the CMB used in this paper, constructed using SEVEN (given at A^side = 64). Other Planck CMB maps used in this 
work are Commader-Ruler, NILC and SMICA, in addition to clean SEVEN maps from 44 to 353 GHz. Right: Planck lensing map, 
optimally filtered to perform the ISW-lensing cross-correlation (given at N^a^ - 1024). 



KSW-bispectrum estimator; bispectra binned in multipole inter- 
vals; and a modal decomposition of the bispectrum. This meas- 
surement is made possible for the first time thanks to the Planck 
data. In addition, we will use information from the lensing field 
as a tracer for an ISW map reconstruction at high redshift (see 
Sect. 6). 

2.2. External data sets 

As described in the introduction, the achromatic nature of the 
ISW eff'ect requires a tracer of the gravitational potentials from 
the large-scale structure, so that by cross-correlating the CMB 
temperature map with that tracer distribution the fluctuations due 
to the ISW effect are singled out. The prerequisites for a tracer 
catalogue to be used in ISW studies are: a large survey volume; 
well-understood biasing properties; and low or at least well- 
modelled systematics. The radio NVSS catalogue and the optical 
Luminous Galaxies (SDSS-CMASS/LOWZ) and Main Galaxy 
Sample (SDSS-MG) catalogues possess these qualities. Table 7 
summarizes some basic properties of these catalogues. In addi- 
tion, the redshift distributions of these catalogues are shown in 
Fig. 2. Notice that NVSS presents the widest redshift coverage. 
The SDSS-CMASS/LOWZ sample is peaked around z ~ 0.5, 
whereas the SDSS-MG sample peaks around z ~ 0.3. 

Figure 3 shows the all-sky density projection for these maps, 
where the grey area indicates regions not observed by these sur- 
veys (or discarded for having contamination or low galaxy num- 
ber density, see next subsections for details). In Fig. 4 we give 
the angular power spectra (blue points) of the surveys (corrected 
with a procedure similar to MASTER, e.g., Hivon et al. 2002), as 
well as the theoretical spectra (black lines) and their 1 <t error 
bars (grey areas), estimated from the MASTER approach as well. 

Besides the cross-correlation between CMB and LSS tracers 
(Sect. 4), we will present results from a different methodology in 
Sect. 5, where we use catalogues of super-structures to study the 
ISW through stacking of the CMB fluctuations on the positions 
of these super-structures. The relevant catalogues are described 
in Sect. 2.2.4. 



2.2.L NVSS radio catalogues 

Luminous Active Galactic Nuclei (hereafter AGN) are known to 
be powerful radio sources, visible out to high redshifts. These 



Table 1. Major characteristics of the galaxy catalogues used 
as tracers of the gravitational potential. From left to right, the 
columns indicate: the number of galaxies per steradian; the frac- 
tion of the sky covered by each survey; the mean bias; and the 
median redshift. Notice that the bias for NVSS is not provided, 
since the assumed model has a bias which depends on redshift 
(see text for details). 



Galaxy catalogue 


n 




/sky 


bias 


z 


NVSS 


1.584 X 




0.73 




1.17 


SDSS-CMASS/LOWZ 


5.558 X 


10'' 


0.22 


2.03 


0.45 


SDSS-MG 


9.680 X 


10* 


0.22 


1.20 


0.32 



N 
"D 

C 
"D 

C 

g 
'4-1 



.12 o 
T3 ^ 



o 
o 





K SDSS-CMASS/LOWZ 




/ \ SDSS-MG 




1 NVSS 





0.0 



0.5 



LO 
redshift z 



L5 



2.0 



Fig. 2. Redshift distributions of the different surveys used in 
this work as LSS tracers, to be correlated with the Planck CMB 
maps. For ease of comparison, these distributions have been nor- 
malised to unity. 



sources are hence able to probe the cosmic density field during 
the entire redshift range from matter domination to accelerated 
expansion due to dark energy. If AGN are fair tracers of the un- 
derlying density field, these sources should likewise probe the 
spatial distribution of the large-scale potential wells that decay 
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-0.5 
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Fig. 3. Density contrast maps obtained from the galaxy cat- 
alogues at A^side = 64. From top to bottom: NVSS; SDSS- 
CMASS/LOWZ; and SDSS-MG. 



at late times after the accelerated expansion sets in and generates 
the ISW effect. 

We shall focus on a single radio survey, with the level of 
sensitivity and sky coverage required for ISW studies, namely 
the NRAO VLA Sky Survey (hereafter NVSS, Condon et al. 
1998). This survey was conducted using the Very Large Array 
(VLA) at 1.4 GHz, and covers up to an equatorial latitude of 

- -40°, with an average noise level of 0.45 mJy beam"'. It 
results in roughly 1.4 x 10^ sources above a flux threshold of 
2.5 mJy. Fig. 3 displays the number density map computed from 
the NVSS survey (top panel). The AGN population is known to 
be dominant in radio catalogues at 1 .4 GHz in the high flux den- 
sity regime. Condon et al. (1998) show that at this frequency, 
star-forming galaxies (SFGs) contribute about 30% of the total 
number of weighted source counts above 1 mJy, but their pres- 
ence decreases rapidly as higher flux thresholds are adopted. The 
NVSS SFGs are nearby sources (z < 0.01), and hence may dis- 
tort the ability of our radio template to probe the intermediate 
and high redshift density field. 
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Fig. 4. Angular power spectra from the maps in Fig. 3. From top 
to bottom: NVSS; SDSS-CMASS/LOWZ; and SDSS-MG. The 
observed spectra are the red points, whereas the theoretical mod- 
els are represented by the black lines (the grey areas correspond 
to the sampling variance). 



We next address the presence of systematic effects in the 
NVSS survey. Two different antenna configurations were used 
while conducting the NVSS survey: the D-configuration (for 
Be e [-10°, 78°]), and the DnC-configuration for large zenith 
angles (b^ < -10°, > 78°). This change in the antenna 
configuration is known to introduce changes in the source num- 
ber density above 2.5 mJy, as first pointed out by Blake & Wall 
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(2002). The NVSS map at 2.5 mJy is corrected for this declina- 
tion systematic using the following procedure: the sky is divided 
into equatorial strips and the mean number of sources in each 
strip is re-normalised to the full sky mean (see e.g., Vielva et al. 
2006). With this procedure the average number of sources in the 
NVSS map is the same as before the correction, and hence the 
shot noise level does not change. The number of strips into which 
the map is divided is 70, but the results are independent of this 
choice. 

Regarding the galaxy bias, in this work we adopt the 
Gaussian bias evolution model of Xia et al. (2011). If n(M,z) 
is the halo mass function and b(M,z) is the bias of halos with 
comoving mass M, then the bias of the survey is given by a mass- 
weighted integral. 



b(z) 



C dM b(M,z) M n(M,z) 
/J dMMn{M,z) 



(2) 



This model depends on the minimum mass M^i„ of halos present 
in the survey. The upper limit in the mass is taken to be infin- 
ity because the effect of the high mass end on the bias is neg- 
ligible. Marcos-Caballero et al. (2013) proposed a theoretical 
model for the NVSS angular power spectrum, which also takes 
into account the information of the redshift distribution given by 
CENSORS data (Brookes et al. 2008). The redshift distribution 
is parametrized by 



dn 
d^ 



"0 



Uo 



(3) 



where zq = 0.32 and a - 0.36. The parameter no is a constant in 
order to have a distribution normalized to unity. This function is 
represented in Fig. 2. The bias follows the prescription of Eq. 2, 
with Mmin equal to 10'^ ''''Mo, where the Sheth-Tormen (Sheth 
& Tormen 1999) mass function is adopted. Hereafter this model 
will be regarded as our fiducial model for NVSS. 

2.2.2. SDSS Luminous Galaxies 

For this analysis we use the photometric Luminous Galaxy (LG) 
catalogue from the Baryonic Oscillation Spectroscopic Survey 
(BOSS) of the SDSS III. The data used consist of two sub- 
samples: CMASS; and LOWZ. Both samples are combined to 
form a unique LG map (see Fig. 3, second panel). Hereafter, 
these samples will be referred to as SDSS-CMASS, SDSS- 
LOWZ, and SDSS-CMASS/LOWZ, for the combination. 

SDSS-CMASS 



We use the BOSS targets chosen to have roughly constant stel- 
lar mass and known as the photometric "CMASS" sample. This 
sample is mostly contained in the redshift range z = 0.4-0.7, 
with a galaxy number density close to 1 10 deg" , and is selected 
after applying the colour cuts explained in Ross et al. (201 1). 

While such color selection yields a catalogue of about 
1,600 , 000 galaxies, farther cuts needed to be applied in order to 
account for dust extinction (based on the maps by Schlegel et al. 
1998 with the criterion E(B - V) < 0.08), for seeing in the r 
band (required to be < 2.0") and for the presence of bright stars, 
similar to Ho et al. (2012). Finally, we neglected all pixels with a 
mask value inferred from the footprint below 0.9 on a HEALPix 
map of resolution A^side=64. This procedure left about one mil- 
lion sources 10,500 deg^. Photometric redshifts of this sample 
are calibrated using a selection of about 100,000 BOSS spectra 



as a training sample for the photometric catalogue. These LGs 
are among the most luminous galaxies in the Universe and there- 
fore allow for a good sampling of the largest scales. Given the 
large number of such sources included in the sample, shot-noise 
does not dominate clustering errors. According to Ross et al. 
(2011), about 3.7% of these objects are either stars or quasars, 
and this makes further corrections necessary, as explained at the 
end of this section. 

SDSS-LOWZ 

The photometric LOWZ sample is one of the two galaxy sam- 
ples targeted by the BOSS of Sloan III. It selectd luminous, 
highly biased, mostly red galaxies, placed at an average red- 
shift of <z> ~ 0.3 and below the redshifts of the CMASS sample 
(z < 0.4). Our selection criteria in terms of the Sloan five model 
magnitudes ugriz follow those given in Sect. 2 of Parejko et al. 
(2013). With a total number of sources close to 600,000, this 
photometric sample contains a higher number density of galax- 
ies in the southern part of the footprint than in the northern one 
(by more than 3 %), which seems to be at odds with ACDM pre- 
dictions. However, most of this effect vanishes when we subtract 
the dipole in the effective area under analysis, in such a way 
that the low £ range of the auto power spectrum is consistent 
with a ACDM model and a constant bias b ^ 2 (Hernandez- 
Monteagudo et al. 2013). 

Both SDSS-CMASS and SDSS-LOWZ samples are further 
corrected for any scaling introduced by possible systematics like 
stars, mask value, seeing, sky emission, airmass and dust ex- 
tinction. Following exactly the same procedure as in Hernandez- 
Monteagudo et al. (2013), we find that both LG samples are con- 
taminated by stars, in the sense that the galaxy number density 
decreases in areas with higher star density, since the latter tend 
to "blind" galaxy detection algorithms. 

2.2.3. Main SDSS Galaxy Sample 

We use a sample of photometrically-selected galaxies from 
the SDSS-DR8 catalogue, which covers a total sky area of 
14,555 deg^ (Aihara et al. 2011). The total number of objects 
labelled as galaxies in this data release is 208 million. From this 
catalogue, and following Cabre et al. (2006), we define a sub- 
sample by selecting only objects within the range 18 < r < 21, 
where this r-band model magnitude corrected for extinction. 
Following Giannantonio et al. (2008), we also restrict our sub- 
sample to objects with redshifts 0.1 < z < 0.9, and with mea- 
sured redsfhit errors such that cr, < 0.5z. We rely on the photo- 
metric redshift estimates of the SDSS photo-z primary galaxy ta- 
ble, which have been obtained through a "kd-tree" nearest neigh- 
bour technique, by fitting the spectroscopic objects observed 
with similar colour and inclination angle. The total number of 
galaxies in our final sample is about 42 million, with redshifts 
distributed around a median value of around 0.35, as shown in 
Fig. 2. To avoid possible errors introduced by singularities in 
the photometric redshifts estimates, instead of using the real ob- 
served redshift distribution in our analysis we resort to the ana- 
lytical function 



dn 



dz r{!2f)z'S^' 



(4) 



which is fitted to the data, with parameters m - 1 .5, j6 = 2.3 and 
Zo - 0.34, which are identical to those found by Giannantonio 
et al. (2012). For the galactic bias we use the value b - 1.2, 
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which is found by Giannantonio et al. (2012) by fitting the 
ACDM prediction to the observed auto-correlation function of 
the galaxies, and we adopt their proposed mask. 



T{() -> T({) ■ 



2.2.4. SDSS, super-structures 

Granett et al. (2008b) produced a sample- of 50 superclusters 
and 50 supervoids identified from the Luminous Red Galaxies 
(LRGs) in the SDSS (sixth data release, DR6, Adelman- 
McCarthy et al. 2008) that covers an area of 7500 deg^ on 
the sky. They used publicly available algorithms, based on the 
Voronoi tessellation, to find 2836 superclusters (using VOBOZ, 
VOronoi BOund Zones, Neyrinck et al. 2005) and 631 super- 
voids (using ZOBOV, ZOnes Bordering On Voidness, Neyrinck 
2008) above a 2 o" significance level (defined as the probabil- 
ity of obtaining, in a uniform Poisson point sample, the same 
density contrasts as those of clusters and voids). The 50 super- 
clusters and 50 supervoids they published in their catalogue cor- 
respond to density contrasts of about 3 cr and 3.3 cr respectively. 
They span a redshift range of 0.4 < z < 0.75, with a median of 
around 0.5, and inhabit a volume of about 5 /i"^ Gpc^. These su- 
perstructures can potentially produce measurable ISW signals, 
as suggested in Granett et al. (2008a,b). For each structure, the 
catalogue provides: the position on the sky of the centre; the 
mean and maximum angular distance between the galaxies in 
the structure and its centre; the physical volume; and three dif- 
ferent measures of the density contrast (calculated from all its 
Voronoi cells, from only its over- or under-dense cells, and from 
only its most over- under-dense cell). For the present paper, we 
concentrate on using the supervoid catalogue by Granett et al. 
(2008b), as they can be compared with two other, more recent 
catalogues of voids. 

The second catalogue of cosmic voids that we consider here 
is published by Pan et al. (2012)''. It has been built from the sev- 
enth data release (DR7) of the SDSS. Using the VoidFinder 
algorithm (Hoyle & Vogeley 2002), they identified 1055 voids 
with redshifts smaller than z - 0.1. Each void is listed with its 
position on the sky, its physical radius (defined as the radius of 
the maximal sphere enclosing the void), an effective radius de- 
fined as the radius of a sphere of the same volume, its physical 
distance to us, its volume and mean density contrast. The fill- 
ing factor of the voids in the sample volume is 62 %. The largest 
void is just over 47 Mpc in effective radius, while the median ef- 
fective radius of the void sample is roughly 25 Mpc. Some of the 
voids are both very close to us and relatively large (larger than 
30 Mpc in radius), which results in large angular sizes of up to 
15°. 

The third void catalogue that we use has been released by 
Sutter et al. (2012) and also made publicly available."^ Note that it 
is being updated regularly, and the results reported here are based 
on the 21 February 2013 version of the catalogue. Using their 
own improved version of ZOBOV, these authors identified 1495 
voids distributed across the 0-0.44 redshift range. They subdi- 
vided their catalogue into six subsamples: diml, dim2, bright] h^c^City = ^ (^"If"?,"'^) 
and bright!, constructed from the main SDSS; and Irgdim, Irg- 
bright built from the SDSS LRG sample. For each void, the in 
formation provided includes the position of the centre, the red 
shift, the volume, the effective radius, and the density contrast. 



3. ISW-lensing bispectrum 

There is an interesting interplay between gravitational lensing of 
the CMB and the ISW effect, which manifests itself as a non- 
Gaussian feature. CMB-lensing can be described by a convo- 
lution of the CMB -temperature map T with the weak lensing 
potential 0, 



(5) 



The CMB lensing can be measured by a direct estimate of the 
CMB bispectrum, because the bispectrum acquires first order 
terms proportional to the product of two power spectra Cj^C^'^, 

where Cj^ is the lensed temperature power spectrum and cj*^ 
is the temperature-potential cross-power spectrum. The poten- 
tial field (f) and the temperature field T are correlated, because 
<p, which deflects the CMB photons, also gives rise to the ISW 
eff^ect in T (Hu & Okamoto 2002; Seljak & Zaldarriaga 1999; 
Verde & Spergel 2002; Giovi et al. 2003). This secondary bis- 
pectrum contains new information about the cosmological red- 
shift, because it is generated mainly at redshifts larger than unity, 
and biases measurements of the primordial bispectrum. The term 
Cj'^ correlates the CMB temperature on small scales with the 
lensing potential on large scales, and causes the bispectrum to 
assume large amplitudes in the squeezed triangles configura- 
tion (see e.g., Goldberg & Spergel 1999; Seljak & Zaldarriaga 
1999; Hu 2000; Giovi et al. 2003; Okamoto & Hu 2003; Giovi 
& Baccigalupi 2005; Lewis & Challinor 2006; Serra & Cooray 
2008; Mangilh & Verde 2009; Hanson et al. 2009, 2010; Smifli 
& Zaldarriaga 201 1; Lewis et al. 201 1). 

Due to the rotational invariance of the sky, the CMB angular 
bispectrum (flf,,„,afj„„flf,„,3) can be factorized as follows: 



(6) 



where = / d" y;!;'(«)y™^(«)F,';''(n) is the Gaunt-integral 

and b[^c^(, is the so-called reduced bispectrum. In the case where 
the bispectral signal on the sky is due to the ISW-lensing effect, 
be^[^i^ - A^'^fejf"^^^^^^, where A^"^ parametrizes the amplitude of 
the effect and 



Llens-ISW 



fA(A + l)- 



2 + 1) + 



' tilcj;cj7 (7) 



H-(5 permutations) >. 



A more general expression for intensity and polarization can be 
found in Lewis et al. (2011). Estimation of the bispectrum then 
yields a measurement of A^*. 

We can also define an alternative rotationally-invariant re- 
duced bispectrum Bf,^,^, as Bc^i^c, = h] i , be^c2h^ where 



(2^1 + 1)(2^2 + 1)(2^3 + 1) / ti h h 
\n 



(8) 



ht tp : //i f a . hawai i . edu/ cosmowave/ 



^ Available at 
supervoids/. 

^ Available at http : //www. physics . drexel . edu/~pan/. 
Available at http : //www. cosmicvoids .net. 



The interest in Bf,^,^, is that it can be directly estimated from the 
observed map using the expression 



Dobs _ ( 



dQr^,(fi)r,,(n)rf3(Q), 



(9) 
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where the filtered maps T({Q) are defined as 



(10) 



By basically combining the single-/" estimates gobs^^iens-isw 
for A^'^ using inverse variance weighting, the ISW-lensing bis- 
pectrum estimator can be written as (see Planck Collaboration 
XXIV 2013, for more details) 



{B 



lens-ISW 



^2jlens-ISW 2jlens-ISW^ ' 

where the inner product is defined by 



(11) 



(12) 



Here, B''" is a linear correction that has zero average but reduces 
the variance of the estimator in the presence of anisotropic noise 
and a mask. Furthermore, Yt^t^e^ = gt^e.j^hj^f^f^Cc^CcCf,, with 
g being a simple permutation factor {g = 6 when all d are equal, 
g - 2 when two { are equal and g - I otherwise). As in all 
expressions in this section, we have implicitly taken the beam 
and noise of the experiment into accout, e.g., C[ should actually 
be b'^Cf + N[ with b( the beam transfer function and N( the noise 
power spectrum. 

In Eq. 1 1 we have also used the fact that, as discussed in de- 
tail in Planck Collaboration XXIV (2013), the full inverse co- 
variance weighting can be replaced by a diagonal covariance 
term, {C^^a)(„ — > a(,„/Cc, without loss of optimality, if the 
masked regions of the map are filled in with a simple diffusive 
inpainting scheme. 

The normalization of the lensing-ISW estimator in the de- 
nominator of Eq. 1 1 can be replaced by (see e.g. Lewis et al. 
2011) 



Eh 



1 + r"^ 
-I ^ ^ '^t 

2£+l 



(13) 



where r( = C J'* 



Ci parameterizes the deviation from the 



Cauchy-Schwarz relation and F( is given in terms of the ISW- 
lensing bispectrum (see for example Lewis et al. 2011). The 
first term in Eq. 13 corresponds to the Fisher errors assuming 
Gaussian a(„. However, contrary to the null hypothesis that is 
assumed, for example, in the primordial bispectra (Gaussianity), 
there is an actual non-Gaussian signal akeady present in the 
ISW-lensing bispectrum. This guarantees a larger variance for 
the estimators than are included in the additional terms present 
in the previous equations. 

An important issue is the impact of the ISW-lensing bispec- 
trum on estimates of the primordial non-Gaussianity. Assuming 
weak levels of non-Gaussianity and considering both the pri- 
mordial bispectrum 5^"™^ ISW-lensing bispectrum 



lens-ISW 

(if if 3 



, one can compute the expected bias A induced in the 



B 

primordial bispectrum using the formula: 



^pnm _ 



(B 



lens-ISW 



(BP™, BP"™) 



(14) 



with the inner product defined in Eq. 12. Predictions of this bias 
on the primordial /nl for Planck resolution can be seen for ex- 
ample in Hanson et al. (2009), Mangilli & Verde (2009), Smith 



& Zaldarriaga (2011), and Lewis et al. (2011). The most im- 
portant bias is introduced to the local shape and, considering 
^max ~ 2000, is expected to be Ai°'="' ~ 7 (Planck Collaboration 
XXIV 2013). 

3.1. ISW-lensing estimators 

There are several implementations of the optimal estimator given 
in Eq. 11. For their detailed description in the context of Planck 
see Planck Collaboration XXIV (2013); Planck Collaboration 
XVII (2013). We have applied four of these implementations 
to Planck data in order to constrain the ISW-lensing bispec- 
trum. Three of them represent a direct bispectrum estimation: the 
KSW estimator (Komatsu et al. 2005; Creminelli et al. 2006), the 
binned bispectrum (Bucher et al. 2010), and the modal decompo- 
sition (Fergusson et al. 2010). The remaining approach is based 
on a previous estimation of the gravitational lensing potential 
field Lewis et al. (2011). These estimators differ in the imple- 
mentation and approximations that are used in order to compute 
the expression given in Eq. 1 1 , the direct computation of which 
is out of reach of current computing facihties. They will be re- 
viewed in the next subsections. 



3.1.1. Lensing potential reconstruction 

The estimator given in Eq. 6 can be written in terms of the lens- 
ing potential amplitude reconstruction ^ as 



A^'f' = S 



1 V rT-^I 

n4^ ' c 



tm 



Tcm 



c 'n! 



(15) 



where 



can be estimated using a quadratic estimator 

(Okamoto & Hu 2003) and Nf is given in terms of the ISW- 
lensing bispectrum (Lewis et al. 201 1). Therefore, this estimator 
quantifes the amount of cross-correlation between the tempera- 
ture map and the reconstruction of the lensing signal, and most 
of the correlation is found at multipoles below 100. 

3.1.2. KSW-estimator 

The KSW bispectrum estimator (Komatsu et al. 2005) for the 
ISW-lensing signal can be written as 



where S can be computed from data as 



tiMi 12^112 ^3^3 



(16) 



(17) 



[(C ^a)(^,„^(C ^a)(^,n^{C ^a)t^mi 

(C 'fl)^„„3j. 



and (F ') is the inverse of the ISW-lensing Fisher matrix F of 
Eq. 13. Details on the implementation of the KSW estimator for 
the ISW-lensing signal can be found in Mangilli et al. (2013). In 
particular, Eq. 16 takes the form 



A — {F )(5 cubic + linear)? 



(18) 



where 5 cubic is the term that extracts the amplitude information 
from the data contained in the bispectrum, while S unear is a zero- 
mean term that reduces estimator variance when the experimen- 
tal setup breaks rotational invariance, i.e., in the presence of sky 
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cut and anisotropic noise. To estimate A^'^ we used the KSW 
estimator with an implementation of the linear term truncated 
at £m-dx as described in Munshi & Heavens (2010) and Planck 
Collaboration XXIV (2013). 

3.1.3. Binned bispectrum 

The binned bispectrum estimator (Bucher et al. 2010) achieves 
the required computational reduction in determining A^* by bin- 
ning Eq. 1 1. In particular, the maximally filtered maps in Eq. 10 
are replaced by 

+t 

where the A, are suitably chosen intervals (bins) of multipole 
values (chosen in such a way as to minimize the variance of the 
quantities to be estimated). These maps are then used in Eq. 9 to 
obtain the binned observed bispectrum, and analogously for Z?''". 
The bispectrum template giens-isw ^j^^ inverse-variance weights 
V are also binned by summing them over all t values in the bin. 
Finally these binned quantities are inserted in the general expres- 
sion for A^'^ (Eq. 11), with the sum over £ replaced by a sum over 
bin indices /. Since most bispectrum shapes change rather slowly 
(with features on the scale of the acoustic peaks, like the power 
spectrum), the binned estimator works very well, increasing the 
variance only slightly, while achieving an enormous computa- 
tional reduction (from about 2000 multipoles in each of the three 
directions to only about 50 bins). 

3.1.4. Modal bispectra 

Modal decomposition of bispectra has been introduced by 
Fergusson et al. (2010) as a way to compute reduced bispectra 
that uses a diagonalization ansatz such that the shape function 
in Fourier space can be separated, which reduces the dimension- 
ality of the integration. At the same time it greatly reduces the 
complexity of estimating bispectra from data. The separation of 
the bispectrum shape function into coefficients qpix) allows the 
derivation of a filtered map Mp(h, x), 

Zqiix)aem 
Y(Jh), (20) 

em 

from the coefficients Of,,, of the temperature map. With that ex- 
pression, one can obtain a mode expansion coefficient j6, 

Pprs = J J x^dx Mp(h, x)Mr(h, x)Ms(h, x). (21) 

With that decomposition, the estimator of the bispectrum as- 
sumes a particularly simple diagonal shape, 

62 

■5 = — A^ 2^ap,;/3prs, (22) 

prs 

where the Oprs are the equivalent coefficients obtained by per- 
forming the modal decomposition of the theoretical bispec- 
trum shape function. The relation between modal bispectra and 
wavelet bispectra is derived by Regan et al. (2013). 



3.2. Results 

The detection of the ISW effect via the non-Gaussian signal in- 
duced by the gravitational lensing potential is summarized in 
Table 2. We provide the estimates of the ISW-lensing ampli- 
tude A^'^, its uncertainty cr^ and the signal-to-noise obtained 
with the different estimator pipelines described in Sec. 3.1. The 
estimators have been applied to the official Planck CMB maps 
made using C-R, NILC, SEVEM, and SMICA (Planck Collaboration 
XII 2013). The quantity cta is obtained from 200 simulations 
representative of the analysed CMB data maps. These Monte 
Carlo simulations (FFP-6, see Planck Collaboration I 2013) ac- 
count for the expected non-Gaussian ISW-lensing signal, ac- 
cording to the Planck best-fit model, and have been passed 
through the different component separation pipelines, as de- 
scribed in Planck Collaboration XII (2013). lensed simulations 
can be found in Planck Collaboration XVII (2013). The mask 
used in the analysis is the combined Galactic and point source 
common mask (U73, Planck Collaboration 1 2013) with sky frac- 
tion /sky = 0.73. 

The KSW and the T(p estimators show similar sensitivity, 
finding, respectively, A^-A = 0.81 + 0.31 and A^-* = 0.70 + 0.28 
from the SMICA CMB map, which corresponds to a significance 
at about the 2.5 cr level. The modal and binned estimators are 
slightly less optimal, but give consistent results, which is con- 
sistent with the imperfect overlap of the modal estimator tem- 
plates with the ISW-lensing signal; the ISW-lensing bispectrum 
has a rapidly oscillating shape in the squeezed limit and both, 
binned and modal estimates, are better suited (and originally 
implemented) to deal with smooth bispectra of the kind pre- 
dicted by primordial inflationary theories. Since the correlation 
coefficient of the binned and modal ISW-lensing templates rel- 
ative to the actual ISW-lensing bispectrum (Eq. 8) is generally 
0.8 < r < 0.9 (to be compared with r = 0.99, achieved by 
both estimators for local, equilateral and orthogonal primordial 
templates, Planck Collaboration XXIV 2013), the corresponding 
estimator's weights are expected to be about 20 % suboptimal, 
consistent with observations. 

The T(p estimator has also been applied to the specific Planck 
lensing baseline, i.e., the MV map, which is a noise-weighted 
combination of the 217 GHz and 143 GHz channel maps, previ- 
ously cleaned from infrared contamination through subtraction 
of the 857 GHz map, taken as a dust template. From this map 
the lensing potential is recovered and then correlated with that 
potential field in order to estimate the amplitude A^"^. The offi- 
cial baseline adopts a more conservative high-pass filtering, such 
that as only multipoles ( > 10 are considered, and the mask 
with /sky - 0.7 is used. In this case, the ISW-lensing estimate is 
0.78 ± 0.32 (a 2.4 cr detection, where the error bars are obtained 
from 1000 simulations), as reported in the first sub-row for T<p 
in Table 2. The full multipole range is considered in the second 
sub-row, obtaining about 7% better sensitivity. 

Notice that, according to all the estimators, the C-R CMB 
map provides lower significance for ISW-lensing, since its res- 
olution is slightly lower than that of the other maps. NILC and 
SMICA exhibit a somewhat larger detection of the ISW signal, 
since they are the least noisy maps. 

In order to explore the agreement among the different es- 
timators, we performed a validation test based on 200 lensed 
simulations processed through the SMICA pipeline. The results 
are summarized in Table 3. For each pair of statistics, we pro- 
vide the difference in amplitudes estimated for the data (AA^"^), 
the dispersion of the difference of amplitudes obtained from the 
simulations (sa), the ratio between this dispersion and the largest 



9 



Planck Collaboration: The ISW effect with Planck 

Table 2. Amplitudes A^'^, errors cta and significance levels of the non-Gaussianity due to the ISW effect, for all component sepa- 
ration algorithms (C-R, NILC, SEVEM, and SMICA) and all the estimators (potential reconstruction, KSW, binned, and modal). For 
the potential reconstruction case, an additional minimum variance (MV) map has been considered (see Planck Collaboration XVII 
2013 for details). 



Estimator C-R NILC SEVEM SMICA MV 

e>\0 0.52 ±0.33 1.5 0.72 ± 0.30 2.4 0.58 ±0.31 1.9 0.68 ± 0.30 2.3 0.78 ± 0.32 2.4 

T'P £> 2 0.52 ± 0.32 1.6 0.75 ± 0.28 2.7 0.62 ± 0.29 2.1 0.70 ± 0.28 2.5 

KSW 0.75 ± 0.32 2.3 0.85 ± 0.32 2.7 0.68 ± 0.32 2.1 0.81 ±0.31 2.6 

binned 0.80 ± 0.40 2.0 1.03 ± 0.37 2.8 0.83 ± 0.39 2.1 0.91 ± 0.37 2.5 

modal 0.68 ± 0.39 1.7 0.93 ± 0.37 2.5 0.60 ± 0.37 1.6 0.77 ± 0.37 2.1 



Table 3. For each pair of estimators we provide the mean differ- 
ence among the amplitudes estimated from the data (AA^**), the 
dispersion of the differences between the amplitudes estimated 
from the simulations (sjC), the ratio of this dispersion to the larger 
of the corresponding sensitivities (rf), and the correlation coeffi- 
cient (p). 









KSW 


binned 


modal 




AA± 
'7 
P 




-0.11 ±0.10 
0.32 
0.95 


-0.21 ±0.21 
0.56 
0.84 


-0.07 ±0.21 
0.56 
0.84 


KSW 


AA± 

n 
p 


sa 




-0.10 ±0.19 
0.52 
0.86 


0.04 ±0.19 
0.51 

0.87 


binned 


AA± 

P 


sa 






0.14 ±0.15 
0.41 
0.92 



of the corresponding sensitivities (j], according to Table 2), and 
the correlation coefficient (p). As can be seen from the Table, 
the agreement among estimators is good and the discrepancies 
are only around 0.5 cr, which is the expected scatter, given the 
correlation between the weights of different estimators discussed 
above. Overall, the bispectrum estimators provide a larger value 
of the amplitude A^*, as compared to the T(fi estimator 

We have also explored the joint estimation of the two bispec- 
tra that are expected to be found in the data: the ISW-lensing; 
and the residual point sources. A detailed description of the 
non-Gaussian signal coming from point sources can be found in 
Planck Collaboration XXIV (2013). The joint analysis of these 
two signals performed with the KSW estimator, and the binned, 
and modal estimators has shown that the ISW-lensing amplitude 
estimation can be considered almost completely independent of 
the non-Gaussian signal induced by the residual sources, and that 
the two bispectra are nearly perfectly uncorrected. 

There is not a unique way of extracting a single signal-to- 
noise value from Table 2. However, all the estimators show evi- 
dence of ISW-lensing at about the 2.5 cr level. 

Finally, we estimate that the bias introduced by the ISW- 
lensing signal on the estimation of the primordial local shape 
bispectrum (Eq. 14) is A^"™ ^ 7, corresponding to the theoret- 
ical expectation, as described in detail in Planck Collaboration 
XXIV (2013). 

4. Cross-Correlation with surveys 

The ISW effect can be probed through several different ap- 
proaches. Among the ones already explored in the literature, the 
classical test is to study the cross-correlation of the CMB tem- 
perature fluctuations with a tracer of the matter distribution, typ- 



ically a galaxy or cluster catalogue. As mentioned in the intro- 
duction, the correlation of the CMB with LSS tracers was first 
proposed by Crittenden & Turok (1996) as a natural way to am- 
plify the ISW signal, otherwise very much subdominant with re- 
spect to the primordial CMB fluctuations. Indeed, this technique 
led to the first reported detection of the ISW effect (Boughn & 
Crittenden 2004). 

Several methods have been proposed in the literature to study 
statistically the cross-correlation of the CMB fluctuations with 
LSS tracers, and, they can be divided into: real space statis- 
tics (e.g., the cross-correlation function, hereinafter CCF); har- 
monic space statistics (e.g., the cross-angular power spectrum, 
hereinafter CAPS); and wavelet space statistics (e.g., the co- 
variance of the Spherical Mexican Hat Wavelet coefficients, or 
SMHWcov from now on). These statistics are equivalent (in the 
sense of the significance of the ISW detection) under ideal condi- 
tions. However, ISW data analysis presents several problematic 
issues (incomplete sky coverage, selection biases in the LSS cat- 
alogues, foreground residuals in the CMB map, etc.). Hence, the 
use of several different statistical approaches provides a more 
robust framework for studying the ISW-LSS cross-correlation, 
since different statistics may have different sensitivity to these 
systematic effects. The individual methods are described in more 
detail in Sect. 4.1. 

Besides the choice of specific statistical tool, the ISW cross- 
correlation can be studied from two different (and complemen- 
tary) perspectives. On the one hand, we can determine the am- 
plitude of the ISW signal, as well as the corresponding signal- 
to-noise ratio, by comparing the observed cross-correlation to 
the expected one. On the other hand, we can postulate a null hy- 
pothesis (i.e., that there is no correlation between the CMB and 
the LSS tracer) and study the probability of obtaining the ob- 
served cross-correlation. Whereas the former answers a question 
regarding the compatibility of the data with the ISW hypothe- 
sis (and provides an estimation of the signal-to-noise associated 
with the observed signal), the latter tells us how incompatible 
the measured signal is with the no-correlation hypothesis, i.e., 
against the presence of dark energy (assuming that the Universe 
is spatially flat). Obviously, both approaches can be extended 
to account for the cross-correlation signal obtained from sev- 
eral surveys at the same time. These two complementary tests 
are described in detail in Sect. 4.2, with the results presented in 
Sect. 4.3. 



4.1. Cross-correlation statistics 

Let us denote the expected cross-correlation of two signals (x 
and y) by , where a stands for a distance measure (e.g., the 
angular distance between two points in the sky, the multipole 
{ of the harmonic transformation, or the wavelet scale R). For 
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simplicity, we assume that the two signals are given in terms of 
a fluctuation field (i.e., with zero mean and dimensionless). 

This cross-correlation could represent either the CCF, the 
CAPS or the SMHWcov. It has to be understood as a vector of 
Qmax components, where amax is the maximum number of con- 
sidered distances. Obviously, when x = y, represents an 
auto-correlation. The specific forms for and Cfw for the dif- 
ferent cross-correlation statistics (CAPS, CCF, and SMHWcov) 
are given below. 

4.1 .1 . Angular cross-power spectrum 

The angular cross-power spectrum (CAPS) is a natural tool for 
studying the cross-correlation of the CMB fluctuations and trac- 
ers of the LSS. Under certain conditions, it provides a statistical 
tool with uncorrelated (full-sky coverage) or nearly uncorrected 
(binned spectrum for incomplete sky coverage) components. 
Even the unbinned CAPS, estimated on incomplete signals, can 
be easily worked out, since the correlations are mostly related to 
the geometry of the mask. This is the case for the CAPS obtained 
through MASTER approach (e.g., Hivon et al. 2002; Hinshaw et al. 
2003). Another approach is to work in the map domain, making 
use of a quadratic maximum likelihood (QML henceforth) esti- 
mator (Tegmark 1997) for the CAPS (Padmanabhan et al. 2005; 
Schiavon et al. 2012). Such approach is optimal, i.e., leads to 
unbiased estimates for the CAPS with minimum error bars. 



Pseudo angular power spectrum 

Let us denote the CAPS between the CMB field T (p) and an 
LSS tracer G (p) map (where p - (0, (p) represents a given pixel) 
as: C^*^ (i.e., = C^*^ for this cross-correlation estimator). In 
the full-sky case, an optimal estimator of the CAPS is given by: 



C 



iTG 



+{ 

m=-e 



(23) 



where ten, and gfm are the spherical harmonic coeflicients of the 
CMB and the LSS maps, respectively. This CAPS can be seen as 
a vector with fmax components, where /"max is the maximum mul- 
tipole considered in the analysis. Here we adopt 3A^side - 1, which 
suffices for ISW analysis, since it is know that most of the ISW 
signal is contained within ^ < 80 (Afshordi 2004; Hernandez- 
Monteagudo 2008). When a mask Yl{p) is applied to the maps, 
it acts as a weight that modifies the underlying harmonic coeffi- 
cients. Now, we have tfm and gf,„, where 



= JJ d (cos 6) dcf>G (6, <f>)n(0,cf>)Y;,„ (6, cf>). 



(24) 



and Y(,„ (6, (p) are the spherical harmonic functions. In these cir- 
cumstances, the estimator in Eq. 23 is not longer optimal, and is 
referred to as pseudo-CAPS. A nearly optimal estimator is given 
by decoupling the masked CAPS (denoted by Cj'^) through the 
masking kernel B (e.g., Xia et al. 201 1): 



where 



2^-H 1 
4-7T 







(25) 



(26) 



with 7p, the cross-angular power spectrum of the T and G masks. 

The estimator in Eq. 25 is nearly optimal because has 
to be understood as the mean value over an ensemble average 
of skies. Let us point out that when more than a single CAPS 
is considered, for instance when one is interested in the cross- 
correlation of the Planck CMB map with more than one LSS 
tracer map, the CAPS estimator can be seen as a single vector 
with Ntm-dx components, with being the number of surveys. 

It can be shown that the element Ca'.ij of the covariance ma- 
trix of the CAPS estimator in Eq. 25 (for the case of a masked 
sky and for surveys) is given by 



C, 



erjj 



(27) 



'■'(It + 1)' 

where 

and (^Mjj^^^^ is the ({,{') element of the inverse matrix of M for 
surveys / and j fixed, such as 



(28) 



M{t,ij = 



2^-H 1 

47T 



Zj"^"^0 J 



(29) 



Here H'^„ is the angular cross-power spectrum of the two joint 
masks, i.e., the masks resulting from the multiplication of the T 
with Gi and Gj, respectively. The quantities Cf are expected or 
fiducial spectra, is the Poissonian noise of the y survey (de- 
convolved by any beam or pixel filter), and 5ij is the Kronecker 
delta. In Eq. 28, the instrumental noise associated with the CMB 
data has been ignored, since the Planck sensitivity is such that 
the noise contribution on the scales of interest is negligible. 
When more than one survey has poor sky coverage, then the 
complexity of the correlations is not well reflected by the previ- 
ous expression. Therefore, in this paper, we will compute Cccjj 
from coherent CMB and LSS Monte Carlo simulations. For each 
simulation, we generate four independent, Gaussian, and white 
realizations (at A^nside = 64), which are afterwards properly cor- 
related using the expected auto- and cross-correlations of the sig- 
nals. Corresponding Poissonian shot noise realizations are added 
to each survey map. The resulting four maps are masked with the 
corresponding masks (i.e., one for the CMB and one for mask for 
each survey). 

The computation of the CAPS in Eq. 25 is extremely fast 
(especially for the resolutions used in the study of the ISW). 
However, as stated above, it is only a nearly optimal estimator of 
the underlying CAPS. Moreover, its departure from optimality is 
largest at the smallest multipoles (largest scales), where the ISW 
signal is more important. 

The QML angular power spectrum 

The QML method for the power spectrum estimation of tem- 
perature CMB anisotropics was introduced by Tegmark (1997) 
and later extended to polarization by Tegmark & de Oliveira- 
Costa (2001). For an application to temperature and polariza- 
tion to WMAP data see Gruppuso et al. (2009) and Paci et al. 
(2013). The same method was employed to measure the cross- 
correlation between the CMB and LSS in Padmanabhan et al. 
(2005), Ho et al. (2008), and Schiavon et al. (2012). The QML 
method is usually stated to be optimal, since it provides unbi- 
ased estimates and the smallest error bars allowed by the Fisher- 
Cramer-Rao inequality. As a drawback, from the computational 
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point of view, the QML is a very expensive approach. Let us de- 
note the QML estimator of the CAPS between the CMB map T 
and an LSS tracer G (at multipole I) by (i.e., ^ = for 
this cross-correlation estimator). 

A detailed description of the algebra of the QML is given 
in Schiavon et al. (2012). We briefly recall here the basics of the 
CAPS estimator, which is given by: 



C'X' 



(30) 



where X represents any of the following spectra: X - {T, TG, G). 
The vector x has 2A^pix elements (with A^pix being the total num- 
ber of pixels allowed by the joint CMB and LSS mask): the first C 
set of A^pix corresponds to the CMB map, and the second one ac- 
counts for the LSS map. The FL9 ^ is the Fisher matrix defined 



7TG X' 



-Tr 



dC 



and the E matrix is given by 



1 dC 



(31) 



(32) 



The object C = S{Cf) -1- N is the global covariance matrix, in- 
cluding the signal S and noise N contributions, with Cf being the 
fiducial theoretical angular power spectrum. The uncertainties in 
the QML estimates are given by the inverse of the Fisher matrix, 
which includes the correlation among different multipoles. The 
error associated with the shot noise of the galaxy surveys is mod- 
elled in the galaxy submatrix of N. 

The results presented on this paper are based on C™, 
whereas Q^'^ is used as a cross-check, applied to a lower res- 
olution version of the maps of A^side = 32. In addition the max- 
imum multipole considered in this case is ^^ax - 2A^side, which 
has been already verified as a conservative limit for the QML. 



4.1.2. Cross-correlation function 

The cross-correlation function (CCF) is a suitable tool for study- 
ing the ISW effect via cross-correlation of the CMB fluctuations 
and tracers of the LSS, and it has been one of the most exten- 
sively used in this context (e.g., Boughn & Crittenden 2002; 
Giannantonio et al. 2008; Xia et al. 2009). On the one hand, 
the signal only appears at very large scales and, therefore, it is 
sufficient to work at resolutions at which the low performance of 
the CCF (in terms of computational time) is not a serious hand- 
icap. On the other hand, neither the CMB nor the LSS data are 
available with full sky coverage and, in some cases, the geome- 
try of the masks is very complicated: the CCF adapts perfectly 
to the effects of partial sky coverage, since it is defined in real 
space. As a drawback, the Poissonian noise of the galaxy tracer 
appears at the smallest angular scales, where the signal-to-noise 
of the ISW effect is higher for this estimator Therefore, a proper 
characterization of the shot noise is especially important for the 
CCF, in order to obtain a good estimation of the uncertainties. 

Let us denote the CCF between the CMB map T and an LSS 
tracer G (at an angular distance of 6*) as C™(6') (i.e., ^" = C™(6i) 
for this cross-correlation statistic). The CCF estimator is defined 
as 



C™(0) 



Ne 



where the sum runs over all pixels with a given angular separa- 
tion. For each angular bin centred around 9, Ng is the number of 
pixel pairs separated by an angle within the bin. Only the pixels 
allowed by the joint CMB and LSS mask are considered. The 
angular bins used in this work are: 9\ e [0, 1]°; 62 e (1,3]°; 
03 e (3, 5]°; ...; and 06i e (119,121]°. The choice of binning 
does not affect the results significantly. 

The covariance of the CCF estimator can be easily derived 
from the one already computed for the CAPS in Eq. 27. It is 
sufficient to know that the CCF can be expressed in terms of the 
CAPS as 



TG 



(«) = J^^C™^ncos0), 

f=0 



An 



(34) 



where P[ (cos 6) are the Legendre polynomials. Hence, it is 
straightforward to prove that the covariance between the 6 and 
ff components of the CCF for the surveys / and j, respectively, 
is given by 

= X Z ^^^^ ^^^T^^f ^) (cos 9') Cw^j-OS) 



An 



4.1 .3. Wavelet covariance 

Wavelets provide an interesting alternative to more traditional 
tools (e.g., CCF or CAPS) for studying the CMB-LSS correla- 
tion. They exploit the fact that the ISW signal is mostly concen- 
trated at scales of a few degrees (e.g., Afshordi 2004). Wavelets 
are ideal kernels to enhance features with a characteristic size, 
since the wavelet analysis at an appropriate scale R amplifies 
those features over the background. Therefore, wavelets could 
provide most of the signal-to-noise of the ISW effect by just 
analysing a narrow range of scales. They were first proposed for 
the ISW detection by Vielva et al. (2006), where the Spherical 
Mexican Hat Wavelet (SMHW, Martinez-Gonzalez et al. 2002) 
was proposed as the filtering kernel. The basic idea of this ap- 
proach is to estimate the covariance of the SMHW coefficients 
(SMHWcov) as a function of the wavelet scale (see e.g., Vielva 
et al. 2006, for details). Other wavelet kernels can be considered, 
such as needlets (Pietrobon et al. 2006a), directional wavelets 
(McEwen et al. 2007), or steerable wavelets (McEwen et al. 
2008). 

Let us denote the SMHWcov between the CMB map T and 
a LSS tracer G (at a wavelet scale R) as (R), i.e., = 
Q™ (R) for this cross-correlation statistic). The SMHWcov esti- 
mator is defined as 



Q. 



TO 



-(Vpix . 



coj,(R)aJc,(R), 



(36) 



where oit (R) and ojc (R) are the SMHW coefficients for 
the CMB and the LSS at scale R, respectively (note that 
wavelet coefficients are forced to have zero mean on the ob- 
served sky). The scales considered in our study are R = 
{60, 90, 120, 150, 200, 250, 300, 350, 400, 500, 600) in arcmin- 
utes. 

As for the CCF, the covariance of the SMHWcov estima- 
tor can be easily derived from the one already computed for the 
CAPS in Eq. 27. It is sufficient to know that the SMHWcov can 
be expressed in terms of the CAPS as 



(33) Q 



TO 



(R) = 2 — 

f=o 



4-K 



oj] (R) , 



(37) 
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where co[ (R) is the SMHW window function at the scale R. 
Hence, it is straightforward to prove that the covariance between 
the R and R' components of the SMWHcov for surveys / and j, 
respectively, is given by: 

C-^'V = Z Z ^ 4 (^') C-^'V- (38) 



4.2. Cross-correlation tests 

For any of the cross-correlation estimators described above, we 
aim two different statistical tests. First, if the observed cross- 
correlation is given by , then, a simple can be proposed to 
estimate the amplitude A, such that A^^ is the closest solution 
to|?': 

(A) = [If - A^f C,>! [I?' - A^] , (39) 

where C^^o is the covariance matrix (of dimension a^ax X «max) 
of the expected cross-correlation i.e., C^n^^ = (^^a^^a'^.y It 
is straightforward to show that the best-fit amplitude A and its 
dispersion are given by 



(40) 



-1/2 



An analogy with Eq. 39 can be defined for the null hypothesis 
case: 



(41) 



where D^n is the covariance of the cross-correlation of the two 
signals, in the absence of an intrinsic dependence, i.e., when 
^a^' = 0. The ISW signal is very weak and, therefore it is a good 
approximation to assume that that D^.y x Cfw. 

For Gaussian statistics, already provides the direct 

probability of the observed cross-correlation under the null 
hypothesis. However, several non-idealities (sky coverage, sys- 
tematics, foregrounds residuals, etc.) forces is to use alternative 
approaches to estimate the probability. One of the most com- 
mon options is to perform the cross-correlation of survey signal 
y with realistic simulations of the CMB, x, and compute a joint 
statistics (s-g-, X^nuw) for each simulation. The probabiltiy value 
associated with the data will come then then be the fraction of 
simulations having a value of x^nuU equal to or larger than the 
one obtained for the data. Both, C^it and D^it can be derived ei- 
ther analytically or numerically (via simulations). 

The latter approach is computationally expensive, but, in 
some cases, could provide a more accurate defence against 
certain systematics, in particular the incomplete sky coverage. 
There are several options to perform such kind of simulations. 
The standard approach is the one mentioned above, i.e., cross- 
correlation of the LSS map with CMB simulations. This is a very 
robust approach, since it is usually hard to reproduce the sys- 
tematics present in the LSS tracers, but incomplete because the 
LSS is fixed, resulting in a lack of randomness. An alternative 
method is to use a jack-knife test, which unfortunately tends to 
underestimate the errors. Finally, one can produce simulations of 
both the CMB and the LSS, assuming perfect knowledge of the 
properties of both signals, in particular of the LSS field (which, 
as mentioned above, is almost never the case). Comprehensive 
discussions of these approaches are given in Cabre et al. (2007) 
and Giannantonio et al. (2008). 



Table 4. Expected significance A/ cta of the CMB-LSS cross- 
correlation. Values obtained from each survey independently, as 
well as jointly, are given for all the estimators (CAPS, CCF, and 
SMHWcov). 





NVSS 


SDSS CMASS/LOWZ 


SDSS MG 


all 


CAPS 


3.0 


1.9 


0.6 


3.2 


CCF 


3.0 


1.9 


0.6 


3.1 


SMMHWcov 


3.0 


1.9 


0.5 


3.1 



4.3. Results 

In this section we present the results obtained from the cross- 
correlation of the galaxy catalogues described in Sect. 2.2 
(NVSS, SDSS-CMASS/LOWZand SDSS-MG) with the four 
Planck CMB maps presented in Sect. 2.1.1 (C-R, NILC, SEVEM, 
and SMICA). All these maps are analysed at a HEALPix resolu- 
tion of A^side - 64. The cross-correlation estimators described in 
the previous section are applied to all cases. This comprehensive 
analysis will help to achieve a robust estimation of the ISW. 

As already mentioned the covariance among all the com- 
ponents of the estimators are obtained from coherent Gaussian 
simulations of the CMB and the three galaxy catalogues. Since 
we are only considering large-scale effects (above about 1°), the 
same set of CMB simulations are equally valid for the four CMB 
maps, since they are nearly identical on such scales(see Planck 
Collaboration XII 2013, for details). We have used 70,000 co- 
herent Monte Carlo simulation sets (as described in Sect. 4.1.1) 
to compute the correlations; this is enough to characterize the 
covariance. 

The expected signal-to-noise ratio for the ISW effect detec- 
tion is summarized in Table 4. Values for all the cross-correlation 
estimators are given. We consider the case of a survey-by-survey 
detection, as well as the joint analysis of all the surveys. A 
signal-to-noise of about 3 cr is expected for the joint analysis, 
which is actually dominated by the NVSS cross-correlation. This 
is expected, since, firstly, NVSS covers a much larger fraction of 
the sky compared to other surveys, and secondly, it extends over 
a redshift interval ideal for the detection of the ISW signal (e.g. 
Afshordi 2004). 

The dififerences among estimators are not significant, indi- 
cating that none of them is clearly optimal compared with the 
others. To explore this agreement further, we have analysed an 
extra set of 1,000 CMB and LSS clustered simulations, and have 
compared, simulation by simulation, the ISW amplitude estima- 
tion derived for each cross-correlation estimator (Cj (9) 
and Q^*^ (R)). In Table 5 we summarize the comparison. We only 
report values for the joint fit to the ISW amplitude for the three 
surveys. Similar results are found survey by survey. For each pair 
of estimators, we provide the mean difference among the ampli- 
tude estimations (AA), the dispersion of these differences (i^), 
the ratio (rj) of this dispersion to the expected sensitivity (i.e., the 
inverse of the signal-to-noise numbers given in the last column 
of Table 4), and the correlation coefficient (p). It is clear that 
the agreement between estimators is very high and that differ- 
ences are, on average, lower than half the statistical uncertainty 
imposed by the sampling variance. 

We have fitted the observed cross-correlations to the ex- 
pected ISW signal (C™, C™ (0), and (R), see Fig. 5), fol- 
lowing Eq. 41, i.e., allowing for a free amplitude of the expected 
signal. Results are summarized in Table 6. Overall, the ISW de- 
tection is at about the 3 cr level and, as expected, it is clearly 
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Table 6. Amplitudes A, errors cr^ and significances Aj cta of the CMB-LSS cross-correlation (survey by survey and all together) 
due to the ISW effect, for all component separation algorithms (C-R, NILC, SEVEN, and SMICA) for the CAPS, CCF, and SMHWcov 
estimators. 



LSS data 


17 


C- 


R 




NILC 




SEVEM 




SMICA 






CAPS 


0.86 + 


0.33 


2.6 


0.91 


± 0.33 


2.8 


on + n 33 


2.7 




2.7 


NVSS 


CCF 


0.80 ± 


0.33 


2.4 


0.84 


±0.33 


2.5 


0.83 ± 0.33 


2.5 


0.84 ±0.33 


2.5 




SMHWcov 


0.89 ± 


0.34 


2.6 


0.93 


±0.34 


2.8 


0.89 ±0.34 


2.6 


0.92 ± 0.34 


2.7 




CAPS 


0.98 + 


0.52 


1.9 


1.09 


±0.52 


2.1 


1 Ofi + 57 


2.0 




2.1 


SDSS-CMASS/LOWZ 


CCF 


0.81 ± 


0.52 


1.6 


0.91 


±0.52 


1.8 


0.89 ±0.52 


1.7 


0.90 ±0.52 


1.7 




SMHWcov 


0.80 ± 


0.53 


1.5 


0.89 


±0.53 


1.9 


0.87 ± 0.53 


1.6 


0.88 ±0.53 


1.7 




CAPS 


1.31 ± 


0.57 


2.3 


1.43 


±0.57 


2.5 


1.35 ±0.57 


2.4 


1.42 ±0.57 


2.5 


SDSS-MG 


CCF 


1.00 ± 


0.57 


1.8 


1.11 


±0.57 


2.0 


1.10± 0.57 


1.9 


1.10 ±0.57 


1.9 




SMHWcov 


1.03 ± 


0.59 


1.8 


1.18 


±0.59 


2.0 


1.15 ±0.59 


2.0 


1.17 ±0.59 


2.0 




CAPS 


0.84 ± 


0.31 


2.7 


0.91 


±0.31 


2.9 


0.88 ±0.31 


2.0 


0.90 ±0.31 


2.9 


all 


CCF 


0.77 ± 


0.31 


2.5 


0.83 


±0.31 


2.7 


0.82 ±0.31 


2.6 


0.82 ±0.31 


2.7 




SMHWcov 


0.86 ± 


0.32 


2.7 


0.92 


±0.32 


2.9 


0.89 ± 0.32 


2.8 


0.91 ± 0.32 


2.9 



Table 5. For each pair of estimators we provide the mean dif- 
ference among the amplitude estimations (AA), the dispersion 
of these differences (s^), the ratio (77) of this dispersion to the 
expected sensitivity (i.e., the inverse of the signal-to-noise num- 
bers given in the last column of Table 4), and the correlation 
coefficient (p). 







CCF 


SMHWcov 




AA ± sa 


-0.01 ±0.12 


0.06 ± 0.07 


CAPS 


1 


0.36 


0.21 




P 


0.93 


0.98 




AA ± Sa 




0.08 ±0.14 


CAPS 


1 




0.42 




P 




0.92 



dominated by the NVSS signal. There are only small differences 
among estimators and CMB maps (as expected from the above 
discussion), indicating that this is a robust result. Notice that 
all the estimated amplitudes are compatible with unity, within 
the error bars (especially for NVSS and SDSS-CMASSA.OWZ). 
This is an additional validation of how CMB and LSS are mod- 
elled. Values of A deviating significantly from unity would in- 
dicate some tension between the observed cross-correlation and 
the model (in particular on the LSS modelling, which is more 
complex). The CAPS-QML, applied to the SEVEM and NVSS 
(i.e., the survey with the highest signal-to-noise), yields a value 
of A = 0.73 ± 0.33, which is compatible with the CAPS, when 
applied to the same A^side - 32 and /"max = 2A^side resolution 
(A — 0.84 + 0.34). Preliminary tests indicate that running the 
CAPS-QML at A^side = 64 resolution could increase the sensitiv- 
ity for detecting the ISW effect with NVSS by ^ 20%. 

Our results indicate a somewhat smaller signal-to-noise with 
respect to some previous analyses on WMAP data, where several 
(and in some case quite similar) surveys were also considered. 
For instance. Ho et al. (2008) and Giannantonio et al. (2012) 
found 3.7 cr and 4.4 cr detections, respectively. Compatibility 
with the former is below 1 cr, whereas there is more tension 
(around 1 .5 cr) with the latter A fraction of around 0.3 cr of these 
differences could be explained in terms of the comological pa- 
rameters adopted to defined the theoretical expectations. In par- 
ticular, the lower values of Hq and Qa found by Planck (Planck 



Collaboration XVI 2013) with respect to WMAP (e.g., Larson 
et al. 201 1), imply a sensitivity for the ISW ^ 10% smaller The 
rest of the differences come either from the LSS side, or from the 
error characterization, which depends on the presence of a corre- 
lated signal between CMB and LSS simulations (see for instance 
Cabre et al. 2007, for a discussion). Survey modelling is another 
important aspect: besides systematic errors associated with the 
galaxy identification and redshift estimation procedures, there 
are complicated aspects, such as the bias characterization. As 
was mentioned akeady, a strong point of our results, is the ex- 
cellent compatibility between the ISW amplitude estimates with 
respect to the expected value. Whereas our estimation deviates 
by about 0.5 cr from the expected value, the Giannantonio et al. 
(2012) result exceeds it by about 1 cr and Ho et al. (2008) is 
around 2 cr above. 

Nevertheless, the value of the ISW effect that we measure 
by means of NVSS (that because of to its large sky coverage, 
redshift range, and density of galaxies is probably the best cur- 
rent catalogue for studying the ISW effect) is significant, and in 
agreement with previously published results using WMAP. 

We have also studied the ISW signal from the point of view 
of its compatibility with the null hypothesis. We have consid- 
ered in this analysis only the NVSS catalogue, since it pro- 
vides the largest detection of the ISW effect and, therefore, is 
the best of the existing surveys for challenging the null hypothe- 
sis. Probability values are summarized in Table 7. As mentioned 
in Sect. 4.2, there is not a unique way of computing the null hy- 
pothesis. Our approach follows Eq. 4 1 , where D^it was computed 
out of 90,000 CMB simulations that have been cross-correlated 
with the LSS data. This matrix is used to compute x^uu ^^om 
the data. This value is then compared to its distribution for the 
null hypothesis, obtained from 1,000 realistic CMB simulations 
(FFP-6) uncorrected with NVSS, which have been processed in 
the same way as the Planck data set. CAPS provides the small- 
est probability value, but the null hypothesis is rejected at about 
10% only; this result is not unexpected, since an expected result 
since the ISW effect is weak. 

The fact that the CAPS statistic provides tighter limits with 
respect to the CCF and SMHWcov could have been anticipated. 
In our implementation, the CAPS explores the maximum an- 
gular range allowed for a given map, whereas the CCF and 
the SMHWcov approaches are only evaluated at certain an- 
gles/scales. This limitation is not an issue in the analysis devoted 
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Fig. 5. Observed and expected cross-correlation signal versus multipole {, for several surveys and different cross-coiTelation es- 
timators. Columns from left to right correspond to: CAPS; CCF; and SMHWcov. Rows from top to bottom represent: NVSS; 
SDSS-CMASS/LOWZ; and SDSS-MG. On each panel we show the expected cross-correlation (black line) and the +1 cr region 
(grey area). Observed cross-correlations for the diiferent CMB maps are provided: C-R and NILC as green and magenta triangles, 
respectively; SEVEM as red circles; and SMICA as blue squares. 



Table 7. Probability values of the CMB-LSS cross-correlation 
for the NVSS survey under the null hypothesis, for the four com- 
ponent separation methods (SMICA, SEVEM, C-R and NILC) and 
for the CAPS, CCF, and SMHWcov estimators. 



LSS data 


17 


C-R 


NILC 


SEVEM 


SMICA 




CAPS 


0.09 


0.10 


0.10 


0.09 


NVSS 


CCF 


0.33 


0.34 


0.40 


0.33 




SMHWcov 


0.20 


0.23 


0.27 


0.19 



to estimating the ISW amplitude, since these angles/scales are 
suitable for detecting the ISW. However, in order to discard the 
null hypothesis, the longer the number of "distances" the better. 

5. Stacking of large-scale structures 

An alternative approach for measuring the ISW effect in Planck 
maps is to look for an ISW signal directly at the positions of 
positive and/or negative peaks in the potential. Since the ex- 



pected (and observed) signal is very weak, for individual struc- 
tures, a stacking technique needs to be applied. Using the WMAP 
data, it has been shown that CMB maps show hot spots and cold 
spots in the direction of superclusters and supervoids, respec- 
tively (Granett et al. 2008a,b, GR08 hereafter), which appear 
to be barely consistent with the predictions of standard ACDM 
(see also Hernandez-Monteagudo & Smith 2012). These struc- 
tures, which are not yet virialized, are evolving while the CMB 
photons travel across them and this should contribute to the 
ISW effect. We apply here the same approach to the different 
Planck maps, using the catalogues of superstructures introduced 
in Sect. 2.2.4, and we test for the robustness of our findings. We 
first discuss our method and the results obtained using the cata- 
logue provided by GR08, and then present the results obtained 
with the other catalogues. 

5.1. Method 

Our analysis is performed on the SMICA CMB map, although 
we have checked that results are compatible for the other three 
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Fig. 6. Stacked regions of Planck maps corresponding to the locations of the superstructures identified by GR08. From left to right 
we show the images resulting from stacking of the 50 superclusters, the 50 supervoids, and the difference of both. The black circles 
superimposed indicate the angular radius at which the signal-to-noise ratio is maximal. See Fig. 7 for the corresponding temperature 
and photometry profiles, as well as their statistical significance. 



Planck maps. We have also used the cleaned frequency maps 
from SEVEN (see Sect. 2.1.1) for some of the tests. We first re- 
move the monopole and dipole of the maps (outside the U73 
mask), and then apply a compact source mask based on the 
Planck Legacy Point Source Catalogue (Planck Collaboration 
Xll 2013) to remove the contamination from individual point 
sources. 

For the purpose of comparison with the results of GR08, we 
smooth the CMB maps with a common Gaussian kernel of 30 ' 
FWHM. We then project them onto patches around each position 
in the supervoid and supercluster catalogue. The GR08 struc- 
tures have a relatively small size on the sky (a few degrees), but 
the other two catalogues considered here contain many larger 
and closer voids, covering larger angular sizes. Thus we work 
with 30° X 30° CMB patches and choose the pixel size to be 
6 ', so that all voids considered are fully enclosed. We then co- 
add (stack) the maps, taking into account the mask used. On the 
stacked images, we calculate both the radial temperature profile 
and the aperture photometry, to characterize the signal around 
density structures. The temperature profile is obtained by com- 
puting the mean of the pixels in rings of fixed width and in- 
creasing angular radius; in practice, it is calculated for 150 radii 
between 0° and 15°, with a width of A0 = 0?1. The photometi-y 
profile is obtained by applying a compensated filter that sub- 
tracts the average temperature of a ring from the average tem- 
perature within the disk whose radius 9 is the inner radius ring, 
and where the outer radius is chosen to be 6 V2, so that the disc 
and ring have the same area. This should enhance fluctuations 
of typical angular size against fluctuations at smaller or larger 
scales. Aperture photometry results are also provide for at 150 
angles, this time between 0° and 15/ V2 =s 10?6. In addition to 
the monopole and dipole, we also removed from the CMB maps 
the contribution of large scale angular modes, namely 1-2- 
10. These modes correspond to angular scales much larger than 
those of the structures under investigation, and for our purposes 
their only effect is to introduce gradients in the stacked images; 
the high-pass filter essentially stops such gradients getting into 
the stacked map (which is equivalent to removing gradients at 
the end). The contribution of the large-scale angular modes has 
no impact on the aperture photometry profiles, and introduces 
only an offset in the temperature profiles (Ilic et al. 2013). 



In order to estimate the significance of the resulting photom- 
etry and temperature profiles, we follow a Monte Carlo approach 
based on stacked CMB images chosen at random positions. In 
detail, we compute the photometry and the temperature profiles 
for 16 000 sets of 50 CMB patches randomly distributed over 
the SDSS area. We then compare the profiles obtained from the 
stacking at the location of the GR08 superstructures to these ran- 
dom profiles, in order to compute their signal-to-noise ratio. 

5.2. ResuWs 

We show in Fig. 6 the stacked images of the 50 supervoids and 
50 superclusters of GR08 in the Planck map. The correspond- 
ing temperature and photometry profiles, along with their sig- 
nificance levels, are shown in Fig. 7. The first thing to say is 
that, although the signatures are fairly weak, the sign of the 
effect certainly seems to be correct. Using the same catalogue 
and the Planck CMB map, we find reasonable agreement with 
GR08. The maximal photometric decrement, -10.8jt/K (essen- 
tially identical with the -11.3;uK found by GR08), induced by 
supervoids is obtained for a preferred scale of about 3?5 (4° in 
GR08) and a signal-to-noise of 3.3 (3.7 cr in GR08), as shown in 
Fig. 7. Superclusters produce a photometric increment of about 
8.5 yuK (slightly above the 7 yuK in GR08), with a significance of 
3.0 cr (compared with 2.6 cr in GR08) at a slightly larger angle 
of 4?7. Finally, the stack of the combined sample (clusters mi- 
nus voids) gives a temperature deviation of 8.7 yuK, with a sig- 
nal strength of 4.0 cr at 4?1, which is consistent with the values 
reported in GR08. The values of statistical significance for our 
aperture aperture photometry results are closely related to those 
for the temperature profiles. Indeed, as shown in the top panel of 
Fig. 7, the temperature profile for the void stack shows a roughly 
2 cr deficit at small angular radii and a roughly 2 cr excess ex- 
tending to large radii. Since the aperture photometry is essen- 
tially an integral of the temperature profile with a compensated 
filter, it picks up enhanced significance because of the shape of 
the temperature profile. 

As noted previously by several authors (e.g., Hernandez- 
Monteagudo & Smith 2012), the amplitude and shape of the pho- 
tometric profile found for voids and clusters is in tension (around 
2 cr) with the values expected from pure ISW within ACDM. 
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Fig. 7. Temperature (top) and photometry (middle and bottom) 
profiles of the stacked CMB patches at the location of the 50 su- 
pervoids and 50 superclusters of GR08. The lower panel shows 
the combined photometry profile (i.e., the average cluster profile 
minus the average void profile). The significance is represented 
by 1, 2, 3, and 4cr level curves (dashed and dotted lines repre- 
sent positive and negative error bars, respectively). These curves 
represent the dispersion of the 16 000 stacks of 50 CMB patches 
chosen at random positions (for illustration, on the top panel, we 
represent in grey 300 of those random profiles). 



However, it is not straightforward to associate this entire signal 
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Fig. 8. Independence of the signal on the SEVEM frequency chan- 
nels. The aperture photometry profiles measured in the stacked 
patches centred on supervoids (top) and superclusters (bottom) 
are virtually identical for all frequencies. 



with a pure ISW effect. As seen in the Fig. 6, many small-scale 
structures - both cold and hot - are present around the region de- 
lineated by the angular radius at which the signal-to-noise ratio 
of the aperture photometry is maximal. This small-scale struc- 
ture contributes to the amplitude of the photometric decrement, 
but at a few tens of juK, which is incompatible with the ACDM 
predictions for the ISW effect. These are rather simply back- 
ground CMB fluctuations, with their lingering presence due to 
the small number (50) of patches which are used to produce the 
stack. 

It is intriguing that the angular sizes of the catalogued su- 
perstructures are smaller than the angular sizes suggested by the 
photometry profiles. This result is more apparent when we re- 
peat the stacking analysis after rescaling each CMB patch by the 
effective radius of the structure it contains. Since the voids and 
super clusters identified by GR08 are roughly the same size, the 
photometric results are similar after rescaling (-lOyuK for voids 
and 7.9 juK for superclusters). However, the deviations have sig- 
nificance levels of 3.3 cr and 2.7 cr for supervoids and superclus- 
ters, respectively, at angular sizes of 1.3 (voids) and 2.6 (clus- 
ters) times the effective radius of the structures. This mismatch 
could be a result of underestimation of the structure extent the 
ZOBOV and VOBOZ algorithms (as already suggested by GR08) or 
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Fig. 9. Stacked images (left) and photometry profiles (right) ob- 
tained for the voids of Sutter et al. (2012). From top to bottom, 
we show results obtained from stacking the largest 231, 589, and 
936 voids. The black circles superimposed indicate the angular 
radius at which the signal is maximal. The statistical dispersion 
is reduced as we stack more voids. However, the amplitude of the 
"cold spot" at about 0.5 times the eff'ective void radius is mostly 
due to the surrounding "hot shell", which is easiest to see in the 
bottom panel. This is further demonstrated by the 3 cr signal de- 
tected using aperture photometry, seen at radii above 1.2 times 
the eff'ective radius (for the 936 void case). Dashed and dotted 
lines in the right-hand panels represent positive and negative er- 
ror bars, respectively, from 1 to 4 cr. 



because larger potential hills and valleys underlie the detected 
superstructures. Since structure in the potential is related to the 
density field through the Poisson equation, gravitational poten- 
tial features are expected to cover larger scales than structures 
in the density field. Nevertheless, the factor of 2.6 for the case 
of superclusters seems large. It is also true that the GR08 super- 
structures were identified in the LRG subsample of the SDSS, 
and LRGs are known to be biased tracers of the matter den- 
sity field (e.g., Tegmark et al. 2006). This bias could help ex- 
plain why structures are larger than the scales identified by the 
VOBOZ algorithm, although the argument would go in the oppo- 
site sense for the voids. Another way of stating this is (Hunt & 
Sarkar 2010) that the relatively large effect decrement found for 
the GR08 voids can be only be attributed to the ISW effect only 
if the underdensities have been significantly underestimated or 
if the LRGs are under-biased. 

It is therefore difficult to be confident that the signal is due 
entirely to the ISW effect. We know, however, that the ISW sig- 
nal generated by superstructures is expected to be achromatic. 



since it generates a fractional perturbation of the CMB temper- 
ature. In order to check if the signal we measure is indeed inde- 
pendent of frequency, we applied the same technique to Planck 
individual SEVEM cleaned frequency maps from 44 to 353 GHz. 
Lower (higher) frequency maps may be contaminated by radio 
(IR) signals coming from our Galaxy and may thus introduce a 
bias in the measurement. Figure 8 shows the photometry profiles 
of supercluster- and supervoid-stacked maps at 44, 70, 100, 143, 
217, and 353 GHz. The flux measured appears to be quite con- 
stant, which supports the idea that the signal is due to the ISW 
effect induced by structures. In the remainder of this section, we 
therefore apply our analysis only to the SMICA CMB map. 

5.3. Discussion and alternative catalogues 

It should be remembered that although the GR08 structures are 
considered to be good tracers of the cosmic matter distribution 
on scales larger than 10 /z ' Mpc, they are also known for their 
sparsity at these redshifts (z ~ 0.4-0.7). This sparsity could lead 
to biased estimates of the properties of the reconstructed voids, 
in particular their sizes and depths could be biased. Moreover, 
some of the structures overlap on the sky, which could lead to 
systematic effects in the stacking analysis. 

We thus turn to other samples, for example the catalogue of 
Pan et al. (2012), introduced in Sect. 2.2.4. The 1054 statisti- 
cally significant voids it contains are larger than 10/;"' Mpc in 
radius and, with redshifts lower than 0.1, they are much closer 
to us than the structures of GR08. Direct stacking gives only a 
weak signal at about the 1 cr level, which is difficult to reconcile 
with the previous results. This may be due to the inclusion of a 
large number of small voids that could dilute the signal. Also, 
unlike the voids of GR08, the voids of Pan et al. (2012) have a 
large scatter in angular sizes on the sky, from about 2° to 20° 
(e.g., llic et al. 2013). In order to try to enhance the signal, we 
repeated the stacking after rescaling the voids to their effective 
sizes. We also subdivided the catalogue into sub-samples based 
on redshift, radius, and/or angular size. However, none of these 
attempts yielded any statistically significant result, in agreement 
with llic etal. (2013). 

Finally, we applied our procedure to the catalogue of voids 
published by Sutter et al. (2012). These cover a rather extended 
range of angular scales (about 2-10°), and so we rescaled all 
the CMB patches by the effective radius of each void. Stacking 
subsample by subsample (diml, dim2, bright!, bright!, Irgdim, 
Irgbright), does not yield any significant signal. Similarly, when 
stacked together, the entire catalogue does not yield an ISW de- 
tection. However, since the ISW signal is expected to be stronger 
for the largest and closest voids (e.g., Flender et al. 2013) we 
tried starting from the largest void and adding them one by one, 
looking for the optimal number of voids, i.e., that for which the 
signal-to-noise ratio is maximal (see Fig. 10). We found that 
stacking 231, 589 or 936 voids gives roughly the same signal-to- 
noise (2.5 cr, 2.0 cr and 2.2 cr, respectively). However, the more 
voids we stack, the smaller the amplitude of the photometry sig- 
nal (see Fig. 9, this being about -2.0//K for 936 voids, -2.1 //K 
for 589 voids and -4.1 fiK for 231 voids, at an angular size of 
about 0.5 times the common rescaled radius. These amplitudes 
are lower than those found with the 50 GR08 voids, although 
above what is expected from numerical simulations (see e.g., 
Hernandez-Monteagudo & Smith 2012, for higher redshift and 
larger voids), but see also Cai et al. (2013). 

The apparent angular size detected (about 0.5 times the effec- 
tive void radius) in the photometry profile is smaller than that for 
the GR08 voids, but in agreement with expectations from simu- 
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Fig. 10. Estimate of the optimal number of patches/voids to 
stack using the catalogue of Sutter et al. (2012). Starting from 
the largest void and adding one CMB patch at a time to the stack, 
we find at each step the minimum of the aperture photome- 
try profile, and we multiply this value by Va^ to find the largest 
signal-to-noise, assuming that the noise scales roughly as 1 / Va^. 
The vertical axis has been normaUzed to the best signal-to-noise, 
obtained for 23 1 voids. 



lations (Cai et al. 2013), and consistent with the sizes detected 
using WMAP data (Ihc et al. 2013). 

The profiles in Fig. 9 show hints (significance about 2 cr) of 
a positive excess below about 0.2 times the effective void radius. 
This may be partly because the ZOBOV algorithm uses galaxies 
as centres of the tessellation, meaning that the centre has to be 
slightly locally overdense. Fig. 9 also shows positive excess for 
larger apertures, partly caused by the large "hot ring" surround- 
ing the cold feature in the stacked images, which raises the mean 
temperature of the stacked image for discs of radii around 0.8- 
1.2 times the void radius. 



6. ISW map recovery 

In recent years, some effort has been invested, not only to ob- 
tain the statistical cross-correlation signal between the CMB and 
LSS data, but also to recover a map of the ISW signal itself 
(Barreiro et al. 2008, 2013; Francis & Peacock 2010; Dupe et al. 
201 1). In particular, assuming the existence of a correlation be- 
tween the CMB and the gravitational potential, it is possible to 
recover a map of the ISW fluctuations using a filtering method, 
given a tracer of the gravitational potential (e.g., the galaxy cat- 
alogues described in Sect. 2.2) and the CMB fluctuations. Given 
the weakness of the signal, the main objective of this section is 
to provide a qualitative image of the ISW fluctuations for visual 
inspection, and an additional consistency test of the validity of 
the assumed fiducial model, by comparing the statistical prop- 
erties of the recovered and expected signals. In addition, this 
secondary anisotropy map could also be used to study the large- 
scale properties of the CMB, and its possible relation to some 
possible large-angle anomalies found in the Planck data (Planck 
Collaboration XXUI 2013). 



6.1. Method 

We have followed the methodology of Bamiro et al. (2008), 
which applies a linear filter to the CMB and to a gravitational 
potential tracer map, in order to reconstruct an ISW map, as- 
suming that the cross- and auto-spectra of the signals are known. 
This technique has been recently applied to reconstruct the ISW 
map from the WMAP data and NVSS galaxy map (Barreiro et al. 
2013). The filter is implemented in harmonic space and the es- 
timated ISW map S[m at each harmonic mode is given by (see 
Barreiro et al. 2008 for details) 



:gem + 



gen 



(42) 



where L{€) corresponds to the Cholesky decomposition of the 
covariance matrix between the considered tracer of the potential 
and the ISW signal, at each multipole, which satisfies C{() - 
L{€)L7{t). Here dcm and g(,„ are the CMB data and the gravi- 
tational potential tracer map, respectively, and C" is the power 
spectrum of the CMB signal without including the ISW effect. 
If full-sky coverage is not available, the covariance matrix is ob- 
tained from the corresponding pseudo-spectra. It can be shown 
that the expected value of the power spectrum for the recon- 
structed signal is given by 



{Cff{\C{0\ + C^iCi) + \Ci{f 

c^(|C(^)l + cfq) 



(43) 



where \C{()\ is the determinant of the tracer-ISW covariance ma- 
trix at each multipole, and C^^ and are the assumed cross- 
spectrum and gravitational potential tracer spectra, respectively. 
Note that the recovered ISW power spectrum will not contain the 
full ISW signal, since it can only account for the part of the ISW 
signal probed by the tracer being considered. It is also worth not- 
ing that in detail the expected cross-correlation depends on the 
assumed model. However, in practice, given the weakness of the 
signal, it would be difficult to distinguish between two mild vari- 
ants of the standard ACDM model. Nevertheless this approach 
still provides a useful consistency check. 

6.2. Results 

We have applied the filter described above to two dififerent cases: 
combining information from the CMB and the NVSS galaxy cat- 
alogue; and applying the filter to the CMB and the recovered 
lensing potential map described in 2.1.2. Results have been ob- 
tained for the four Planck maps, C-R, NILC, SEVEM, and SMICA. 
For simplicity, we show the reconstructions only for the SEVEM 
CMB map, since the four methods give very similar results. The 
resolution considered for both analyses is A^si^e = 64. 

For the first case, we are using the Planck fiducial model 
for the CMB and cross-power spectrum, while for the NVSS 
map we assume the model described in Sect. 2.2. 1. We also take 
into account the presence of Poissonian noise. We have excluded 
the area obtained from combining the CMB mask at A^side = 64 
(described in Sect. 2.1.1) as well as the area which has not been 
observed by NVSS. The final mask keeps around 62% of the sky. 
Since the filter is constructed in harmonic space, we have used 
an apodized version of the mask in order to reduce the mask- 
induced correlations. In any case, the degradation introduced by 
the presence of a mask is small (Barreiro et al. 2008). 

For the second case, the lensing map involved applying a 
high-pass filter, which removed all multipoles with ^ < 10. This 
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Fig. 11. Reconstructed ISW map from the Planck CMB and NVSS data (left) and from the Planck CMB and lensing potential maps 
(right). Note that the maps are not expected to look exactly the same, since each of them provides a partial reconstruction of the 
noisy ISW signal (see Sect. 6.2 for details). 



filtering was done in harmonic space with the presence of a 
mask. To take this effect into account we used a direct estima- 
tion of the pseudo-power spectrum of these data for the power 
spectrum of the lensing map, after applying the corresponding 
apodized mask. We used the Planck fiducial model for the other 
power spectra involved, but setting to zero the cross-power for 
{ < 10. A mask has been constructed by combining the CMB 
mask plus that provided for the lensing potential map (described 
in Planck Collaboration XVII 2013), which keeps around 67% 
of the sky. The corresponding apodized version of this mask was 
applied before reconstructing the ISW map. Note that the map 
given in Fig. 1 (right panel) corresponds, to a good approxima- 
tion, to the first term of the right hand side of Eq. 42. 

Figure 1 1 shows the reconstructed ISW map using the 
Planck CMB map and NVSS (left panel) and that obtained com- 
bining the CMB with the lensing potential map (right panel). 
There are similar structures present in both maps, but they are 
not expected to look exactly the same, since each of them pro- 
vides only a partial reconstruction of the ISW signal. This is 
due to the fact that the reconstruction accounts for the part of 
the ISW effect probed by the considered tracer, which is dif- 
ferent (although correlated) for each case. Moreover, due to the 
high-pass filter applied to the lensing potential map, the power 
2A. ( < 10 for this case corresponds to the Wiener- filtered map 
of the CMB (to which the filter given by Eq. 42 defaults, if the 
cross-correlation is set to zero, as in this case), without additional 
information from the considered tracer. 

For both cases, we have tested that the power spectrum of the 
recovered ISW signal, as well as that of the cross-power between 
the reconstructed ISW and the considered gravitational potential 
tracer, are consistent with the corresponding expected values. 
This indicates the compatibility between the assumed fiducial 
model and the underlying statistical properties of the data. 



7. Conclusions 

This paper presents the first study of the ISW effect using Planck 
data. We derived results based on three different approaches: the 
detection of the interplay between weak lensing of the CMB and 
the ISW effect, by looking at non-Gaussian signatures; the con- 
ventional cross-correlations with tracers of large-scale structure; 
and aperture photometry on stacks of the CMB field at the po- 
sitions of known superstructures. A reconstruction of the ISW 
map inferred from the CMB and LSS tracers was also provided. 



The correlation with lensing allows, for the first time, the de- 
tection of the ISW effect using only CMB data. This is an effec- 
tive approach, because the gravitational potential responsible for 
deflecting CMB photons also generates ISW temperature per- 
turbations. Using different estimators, we investigated the cor- 
relation of the Planck temperature map with a reconstruction of 
the lensing potential on the one hand, and the estimation of the 
ISW-lensing generated non-Gaussian signature on the other. We 
found that the signal strength is close to 2.5 cr, for several com- 
binations of estimator implementation and foreground-cleaned 
CMB maps. 

We computed cross-correlations between the Planck 
CMB temperature map, and tracers of large-scale structure, 
namely: the NVSS survey of radio sources; and the SDSS- 
CMASS/LOWZ, and SDSS-MG galaxy samples. As estimators 
we considered the angular cross-correlation function, the angular 
cross-spectra, and the variance of wavelet coefficients as a func- 
tion of angular scale. We performed a comparison on different 
component-separation maps, where we considered C-R, NILC, 
SEVEM, and SMICA, and found remarkable agreement between 
the results, indicating that the the low multipoles are robustly re- 
constructed. Covariance matrices between the cross-correlation 
quantities were estimated for a set of Gaussian realizations of the 
CMB for the Planck fiducial model. For the ISW effect, we re- 
port detection significance levels of 2.9 cr (NVSS), 1 .7 cr (SDSS- 
CMASS/LOWZ), and 2.0 cr (SDSS-MG), which are consistent 
among the different estimators considered. Although these num- 
bers are compatible with previous claims which used WMAP 
data, they are generally smaller. We believe that this discrep- 
ancy is mainly due to the different characterization of the sur- 
veys and treatment of uncertainties, since the measurement of 
the CMB fluctuations at the scales which contribute to the ISW 
detection are very similar for Planck and WMAP. Only a fraction 
of these differences (around 0.3 cr) could be understood in terms 
of the different cosmological models used by each experiment - 
in particular, the lower values of Hq and Q./^ reported by Planck 
compared with WMAP. 

A strength of our new study lies in the fact that the ampli- 
tudes derived for the expected signals are largely consistent with 
unity (i.e., the model expectation), which indicates good mod- 
elling of the surveys. The CMB and LSS cross-correlation has 
also been tested against the null hypothesis, i.e., whether the ob- 
served signal is compatible with a null correlation. As expected 
for such a weak signal, there is no strong evidence of incom- 
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patibility with the lack of correlation. In this respect, the CAPS 
approach seems to provide better constraints than the other esti- 
mators investigated here (CCF and the SMHWcov). 

We explored the aperture photometry of stacked CMB 
patches at the positions of superstructures identified in the SDSS 
galaxy distribution. Our analysis of the Granett et al. (2008a) cat- 
alogue (50 supervoids and 50 superstructures) reproduced pre- 
vious results, with similarly strong amplitude and significance 
levels (somewhat above and below 3 cr for voids and clusters, 
respectively). While the most plausible source of this signal is 
the ISW effect associated with these structures, it shows some 
tension with expectations, both in terms of amplitude and scale. 
The same type of analysis was carried on the latest and much 
larger void catalogues of Sutter et al. (2012) (about 1 500 voids) 
and Pan et al. (2012) (about 1 000 voids). The results range from 
negligible to evidence at the 2-2.5 cr level, with a more mod- 
erate amplitude and a smaller scale, in better agreement with 
theoretical predictions found in the literature. The broad spectral 
coverage of Planck allows us to confirm the achromatic nature 
of these signals over the 44 to 353 GHz range, supporting their 
cosmological origin. 

We reconstructed maps of the ISW effect using a linear fil- 
ter, by combining the Planck CMB and a gravitational potential 
tracers. In particular, we considered both the NVSS catalogue 
and the reconstructed CMB lensing map as LSS tracers. Again 
we found good agreement between different component separa- 
tion methods, as well as consistency between the expected and 
reconstructed auto- and cross-power spectra for the recovered 
ISW map. 

We conclude that the ISW effect is present in Planck data 
at the level expected for ACDM-cosmologies, using a range of 
measurement methods, although there is a possible tension with 
the results from stacking of CMB fields centred on superstruc- 
tures. Generally, our results are more conservative than previ- 
ous claims using WMAP data, but the agreement with the ex- 
pected signal is better Future Planck data releases, including 
polarization information, as well as improved understanding of 
foregrounds, could improve on these results, in particular for 
ISW-lensing correlation and ISW-lensing map reconstruction. 
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